Quantum Superposition States of Bose-Einstein condensates 



J. I. Cirac 1 , M. Lewenstein 2 , K. M0lmer 3 and P. Zoller 1 

(1) Institut für Theoretische Physik, Universitat Innsbruck 
Technikerstrasse 25, A-6020 Innsbruck, Àustria 
(2) Commissariat a l'Enegie Atomique, DSM/DRECAM/SPAM 
Centre d'Etudes de Saclay, 91191 Gif-sur-Yvette, France 
(3) Institute of Physics and Astronomy, University of Aarhus, DK-8000 Àrhus C 



Abstract 

We propose a scheme to create a macroscopic "Schòdinger cat" state formed 
by two interacting Bose condensates. In analogy with quantum òptics, where 
the control and engineering of quantum states can be maintained to a large 
extend, we consider the present scheme to be an example of quantum atom 
òptics at work. 
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I. INTRODUCTION 



The recent experimental realization of Bose-Einstein condensation of trapped cold ru- 
bidium U, sodium |J, and lithium || àtoms has initiated new areas of atòmic, molecular 
and optical physics ||. While some of these new areas remain still somewhat speculative, 
others have already attained firm experimental grounds, and many of them are based on the 
analogy between the matter waves and electromagnètic waves, or between bosonic àtoms 
and photons. 

On the level of single àtoms, the analogy between the matter and electromagnètic waves 
has led to ràpid developments of atom òptics ||. Some authors have thus considered a 
possibility of nonlinear atom òptics in systems of many cold àtoms, where the quantum 
statistical properties and atom-atom interactions become important ||. It has been also 
pointed out [0] that nonlinear excitations of Bose-Einstein condensates (BEC) may lead to 
various analogs of nonlinear òptics. Most of these theories have a mean field character, i.e. 
they disregard quantum fluctuations of the atòmic field operators and concentrate on the 
nonlinear Schòdinger-like wave equations for the matter wave functions. 

In many situations, such as for example làser cooling or optical imaging, cold àtoms 
not only exhibit their quantum statistical properties, but on top of that they interact with 
photons. This fact motivated the developments of quantum field theory of àtoms and photons 
||. Although this theory accounts in principle for quantum fluctuations of both atòmic and 
electromagnètic fields, the attention has so far been focused predominantly on the latter. 

The atom-photon analogy has also triggered the studies in the area of physics of atom 
làsers |§. An atom làser, or a boser is a matter wave analog of an ordinary làser. Quite 
recently, the possibility of employing BEC as a source of coherent matter waves has been 



demonstrated in the remar kable experiments of the MIT group [10|. 

We propose here to proceed with this analogy and to look for the implementation of 
further elements of quantum òptics in quantum atom òptics. In our opinion, one of the 
major domains of concern of modern quantum òptics is preparation, engineering, control and 



detection of quantum states in various systems that involve light-matter interactions |l I 
By analogy, quantum atom òptics, in the sense proposed, concerns preparation, engineering, 
control and detection of quantum states in atòmic systems. Recent studies of excitations in 
trapped BEC belong to this category, although so far rather limited kinds of time dependent 
perturbations of the trapping potential ]T2] were used and only few types of excitation have 
been investigated. Walsworth and You have proposed a method of selective creation of 
quasi-particle excitations in trapped BEC [[OJ . Their method, referred to as spatial magnètic 
resonance method, could in principle allow for engineering and control of arbitrary excitation 
in the Bose condensed system. 

One of the most spectacular achievements of quantum òptics in the recent years has been 
the observation and study of macroscopic (or strictly speaking mesoscopic) "Schròdinger 



cat" states of a trapped ion |14|, and of an electromagnètic field in a high Q cavity ||Ï5 
Schòdinger |TB[ has introduced his famous "cat" states in order to illustrate the fundamental 
problem of the correspondence between the micro- and macro-worlds: the fact that quantum 
superposition states are never observed on the macro level. As postulated by von Neumann 
|T7|j , this is due to irreversible reduction of superposition states into statistical mixtures. Such 
reduction occurs in any quantum measurement process and leaves the considered system in 
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a mixed state in a "preferred" basis, determined by the measurement. Modern quantum 
measurement theory [18| describes the reduction process in terms of quantum decoherence 
due to interactions of the system with environment | ÏÜ| . Experimental realization of "cat" 



states requires thus typically sophisticated means to avoid the decoherence effects []T4| , |T5 



In this paper we demonstrate that it is feasible to prepare, control and detect a 
"Schròdinger cat" formed by two interacting Bose condensates of àtoms in different internal 
states. Atom-atom interactions in our model are mediated through atom-atom col·lisions and 
a Josephson-like làser coupling that interchanges internal atòmic states in a coherent man- 
ner. The theory of such bi-condensates restricted to collisional interactions only has been 



recently discussed in the Thomas- Fermi approximation ||20|| , and beyond pl|j2^] . Amaz- 
ingly, the simultaneous condensation of 87 Rb àtoms in two internal states (F=2,M=2) and 
(F=1,M=-1) has been recently achieved at JILA ||23|| , using a combination of evaporative 
HH and sympathetic |25| cooling. As pointed out by Julienne et al. f2(|, the simultaneous 
condensation was possible due to a very fortunate ratio of elastic/inelastic collision rates for 
Rb àtoms. For the moment, the perspectives of extending the result of |23[ to other atòmic 
species are not promising. Nevertheless, one could expect that various ways of modifying 
atòmic scattering lengths will be realized [^7]], which will allow to control atòmic collision 
processes in a desired way The above comments apply to the case of magnètic traps. Once 
it becomes possible to achieve Bose-Einstein condenstaion in, for example, far-off resonance 
traps, this will open other possibilities of trapping partides with different internal levels. 
This will allow one to meet in real atòmic systems the conditions discussed below for the 
preparation of Schròdinger cat states. 

The plan of the paper is the following. In section II we present the quantum field 
theoretical model of two trapped condensates, and its simplified two-mode caricature. The 
detailed analysis of the two-mode model is carried out in section III. We show that in 
some circumstances the ground state of the system becomes a "Schròdinger cat" , and that 
the system can be prepared in such a state by adiabatically changing the strength of the 
Josephson-like làser coupling. Various approximate solutions are tested here in comparison 
with the exact numerical solution of the problem. In section IV approximate solutions of 
the complete quantum field theoretical model are found. They display the same qualitative 
behavior as the one obtained for the two mode model. Finally, in section V we discuss the 
feasibility of experimental observation of "Schròdinger cat" states of two interacting Bose 
condensates. 



II. QUANTUM FIELD THEORY OF TWO INTERACTING CONDENSATES 

We consider here Bose-Einstein condensation of a trapped gas of àtoms with two internal 
levels \A) and \B). The àtoms interact via AA, BB, and AB elàstic col·lisions. Additionally, 
a set of làser fields drives coherently a Raman transition connecting \A) <-> \B). In the 
formalism of second quantization, such a system is described by the following Hamiltonian 

H = H A + H B + Hi nt + H Las (1) 

where 
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Here, Ha, b describes the evolution of àtoms in \A) and \B), respectively, in the absence 
of interactions between àtoms in different internal states; iïi nt describes the interactions 
between àtoms in \A) and \B) due to col·lisions; i?Las describes the Raman transitions induced 
by the làser detuned by A from the Raman resonance; such interactions act as a Josephson- 
like coupling which transfers coherently partides between \A) and \B), at a Rabi frequency 
O > 0. Àtoms are confined in harmònic potentials of freqüències u> a ,b, and 

q>a b ab are 

the scattering lengths for the corresponding col·lisions, respectively We assume that the 
col·lisions are purely elàstic, and that they do not change the number of partides in each 
internal level. 

The field operators ^a,b{x), í&a,b(%) annihilate and create àtoms at x in the internal 
state \A) and \B). They fulfill the Standard bosonic commutation relations 



[^a{x)^a{x') ] ] = 8{x-x' 



[*b(^),*b(5 
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For the sake of simplicity, throughout this paper we will assume that 
làser excitation A = 0, and u B = uja = oj. This makes the Hamiltonian (0) invariant 
under the exchange A B, which simplifies the analytical arguments. In experiment 
(c.f. ||23|| ) this symmetric situation is not directly realized, since àtoms in different (F,M) 
states experience different Zeeman effects in the magnètic field, and feel thus different trap 
potentials. Nevertheless, we stress that our assumption has only a technical character. 
All the results presented below for the A — B symmetric case can be translated into the 
asymmetric case, as we shall see below. In general, if uja ^ u B , one can always choose the 
detuning A ^ to compensate for the potential difference. One should also mention the 
fact that the different Zeeman effects, combined with gravity, may displace the traps with 
respect to each other; even for such a sitation compensation of potential difference using 
A ^ should be posible, although more complex. 

Let us now rescale the variables to dimensionless ones as follows: First, we divide H 
by hu), and then define A = Q/(2uj) > 0, f = x/x , 4f(r) = &(x)x^ 2 , U = ATca sc /x , and 
Ui = 4na s Jl B /xo, where xq = (h/Mcu)^ 1 ^ 2 ^. The rescaled dimensionless Hamiltonians (0) read 
now 



H a ,b = J d ó r^A, B (r 
H lnt = U x j d 3 f^ A (r)^B(r)^B(r)^A 
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Our objective is to study the properties of the system described by the above Hamilto- 
nians at zero temperature. To tiris aim we need to find the ground state of Hamiltonian 
(0) with the corresponding terms defined in Eqs. (H). The search of this state is a very 
difficult task. Under some conditions one can obtain mean field approximations, and nu- 
merical (approximate) results for the ground state of ([!]). In order to understand better our 
model, we will first analyze a very simple two-mode model described by a caricature of the 
Hamiltonian (|Ip. As we shall see in section III and IV, the analysis of the simplified model 
resembles very much the analysis of the complete model described by ([[]). 

The two-mode approximation of the Hamiltonian ([TJ) is given by (|T|) but with 

Ha = oja^ol + —^-afafaa, (5a) 

H B = u) B b ] b + ^tftfbb, (5b) 

Hint = Uiaïtfba, (5c) 
#Las = -A (<Jbe~ iM + tfae +iAt ) . (5d) 

This model corresponds to the previous one in the limit in which the external motion of the 
àtoms is frozen. The bosonic annihilation and creation operators a,a^ , b, b' annihilate and 
create àtoms in infernal states \A), and \B), respectively. They fulfill Standard commutation 
relations= [a,a^\ = [b,b'] = 1. As before, we will assume ujb = = A = 0, and 
U aa = Ubb = Uq. This allows us to neglect the first term in Ha and H B , since the total 
number of partides N = a^a + tfb is conserved. Note that the same simplification occurs 
when uja ^ ^b, but A = uja — oüb- This means that the results obtained in the next 
section for A — B symmetric CclSG, £1X6 equivalent to the ones for the asymmetric case with 
appropriately chosen A. 

An additional motivation behind the use of the model @ is that it can be solved nu- 
merically for moderate N, and therefore it allows to compare the analytical approximations 
with the exact numerical results. This will provide us with a guideline for the analysis of 
the complete quantum field theoretical model section III. 



III. ANALYSIS OF THE TWO-MODE MODEL 

In this section we consider in details the simple two-mode model described by the Hamil- 
tonians (|5|). The section is divided into 3 subsections. In the first of them we derive the 
ground state energy of @ using a mean field approach, for which all the àtoms are supposed 
to be in the same single particle state. In subsection III. B we refine this theory to find a bet- 
ter approximation to the ground state. We show that under certain conditions the ground 
state corresponds to a Schròdínger cat state. Finally, in subsection III. C we diagonalize the 
Hamiltonian exactly using a numerical method (|5]), and compare the exact predictions with 
the approximate ones of the previous subsections. 
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A. Mean Field Approximation 



The equations for the ground state in the mean field (Hartree-like) approximation can 
be derived using the Standard procedure. We consider the single particle state 



\il> 1 ) = a 1 \A)+p 1 \B), 



(6) 



where |ai| 2 + \/3i\ 2 = 1, and look for collective states of N particle system, with all the par- 
tides in the same state (0) which minimizes the total energy. Using the second quantization 
description, these collective states can be represented as 



\i/>n) = l^i) ® \i/>i)...\il)i) 



1 



N 



+ (3 x tf |0) 



(7) 



where |0) denotes the vacuum state. 

The expectation value of the Hamiltonian (El) in this state is 



E(a, a*, (3, (3*) = (^ N \H\^ N ) = ^ (|«| 4 + |/?| 4 



C/l I Q! I 



where we have redefined £7 ,i = U 0j i(N — 1)/N, and a 
normalization condition imposes now 



- A {a*f3 + f3*a) , 
A^ai and (3 = y/Nfo. 



\a\ 2 + 



N. 



(8) 
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For simplicity of the notation, we will drop the tilde over U's in the following. 

We minimize the mean energy (§) with respect to a,/3 and their complex conjugates, 
imposing the constraint @ by using a Lagrange multiplier fi. After elementary calculations 
we obtain, 



U \a\ 2 + U^Pl 2 a — \p = ya, 
U \f3\ 2 + U^al 2 } (3 - Xa = /2(3. 



(10a) 
(10b) 



The above equations can be easily solved. To this aim, we first note that for A ^ 0, a and 
(3 can be taken to be nonvanishing real numbers. Thus, we can divide the first equation by 
a and the second by (3, and subtract them to obtain 



a/3 



(I«| 2 -|/3P 



0. 



(11) 



The analysis of Eq. ( |ïl~|) is straightforward. Defining A = 2\/\N(U\ — Uq)\ one finds that 
for |A| > 1 there exists only one solution 



«o = (3 = 



which gives the mean energy energy 

E = 



N 2 



XN: 



(12) 



(13) 
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For |A| < 1 there exist three solutions (0, +, — ): 



a± = P T 
with the corresponding energies 

E = 



a =(3 = v /ÏV/2 
■N 



(1±VT 



1/2 



AT2 

Ei = E_ 



(U + C/i) - AA^ 
TT N 2 X 2 

— Un — 



Ux-Uo 



(14a) 
(14b) 

(15a) 
(15b) 



One can easily check that for Un > U\ we have E± > E and therefore the solution (an, fio) 
gives the minimum energy. On the other hand, for U\ > Un both solutions (a + ,/3 + ) and 
(a_,/3_) give a lower energy than (a , (3 ) (in particular, for A = 1, E — E±). 
The results can be summarized as follows: 

(a) Weak interactions between àtoms in the states \A) and \B): In this situation Un > U\ 
and the mean field wavefunction for the ground state is 



1 



N 



+ 6 f |0) 



Vn n n\ 



N 



anaï + pntf |0). 



(16) 



(b) Strong interactions between àtoms in the states \A) and \B): In this situation U\ > Un 
and one has to distinguish two cases: 

• (b.l) Strong làser case: A > 1 and the mean field wavefunction is given in 
(0)- 

• (b.2) Weak làser case: A < 1; there are two degenerate solutions |^) for the 
mean field ground state wavefunction 



a±a ] + (3±tf |0) 



Vn n n\ 

where a± and (3± are given by expression (|14b|) 



(17) 



B. Beyond the Mean Field Approximation 



For the chosen parameters, the Hamiltonian ([!]), with (^|) is invariant under the operation 
Tab that exchanges the infernal level \A) with \B). Thus, in case of no degeneracy the 
eigenstates of H, \<f>k) must be eigenstates of Tab too. Since Tab is idempotent (i.e. has 
eigenvalues ±1), the eigenstates have to fulfill Tab^u) = ilV'fe)- For \a\ ^ \(3\ (the case b.2 
above), it is clear that the states obtained using the mean field approach (|7|) do not satisfy 
this condition. This indicates that in the case (b.2) one can obtain a better approximation 
to the ground state with a lower energy if one uses as a variational ansatz the wavefunction 



7 



m = (ii&> ± m) /V2. (is) 

This is a superposition of the two degenerate solutions. Note that |V' ± ) is indeed an eigenstate 
of Tab with eigenvalue ±1, and therefore it conforms with the symmetry of the Hamiltonian. 

The states (|I8|) are written as a superposition of two states in which either all the àtoms 
are in the single particle state j^i ) = aii\A) + (3i\B), or all are in the single particle 
state \ipï) = aï\A) + (3ï\B). Therefore, they have the form of Schròdinger cat states. 
Note, however, that a Schròdinger cat state is characterized by its coherent inclusion of 
macroscopically distinguishable states. For the state of our condensates to be a true (i.e. 
macroscopic, or at least mezoscopic) Schròdinger cat state we must therefore require that 
the overlap e, 

e = (^\^)=A N , (19) 

be as small as possible. The "size of the cat", which can be defined as 1/e, should, on the 
other hand be as large as possible. The theory should determine under which conditions the 
observation of the "cat of maximal size" is feasible. 

The expectation value of the energy of the state is 



N 2 2U - A 2 (f/i - U ) ± (A) N (3U - U x ) 



(20) 



4 1 ± (A) N 

It is easy to check that in the limit of e ~ 1 (i.e. when the cat is still microscopic) , we obtain 

N 

AE = -(U 1 -U ). (21) 

This Eq. reveals characteristic scaling of the energy difference AE with N, which, as we 
shall see below, is also vàlid in the more interesting limit of e = <C 1 (i.e; when the cat 
is mezo-, or macroscopic). In such case we may first expand the result ( p0|) in e, and then 
in A" 3> 1, so that we obtain 

AE = E_-E + -~e\n(e)N(U 1 -U ). (22) 

Thus, for a "given size" of the cat 1/e the energy difference is proportional to N. Quite 
generally, the difhculty in cooling to a ground state of a given purity increases with growing 
number of àtoms N, while a larger energy gap AE makes the cooling easier. In this sense, 
the scaling AE oc A^ helps. 



C. Exact Numerical Solution 

In the Fock basis \m)A ® \N — m)B (ui = 0, 1, ... , N) the Hamiltonian H is a (N + 
1) x (N + 1) tridiagonal matrix, and can therefore be easily diagonalized by numerical 
methods. Since the mean field approximation and its improved version analyzed in previous 
subsections should be vàlid in the limit A^ — > oo, we concentrate here on the finite A^ results. 
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Let us denote by 

N 

\<Pk)= Y. ( t\^)A®\N-m) B , (23) 

m=0 

the eigenstate corresponding to the energy Ef. (k — 0, 1, . . . , N, and E < E\ < . . . < E^). 
The results of our analysis are presented in Figs. 1-4. In Fig. 1 we have plotted the ground 
state energy Eq as a function of A for N = 1000 and U\ = 3Uq. Although this Figure 
already shows the clear signatures of the "phase transition" to the "Schròdinger cat" phase 
for A < 1, it is more instructive to look at relative behavior of the consecutive eigenenergies 
of the low excited states. This is represented in Fig. 2(a) for N = 1000 and 2(b) for 
N = 10000, where we have plotted the ratio between the energy difference of the first 
excited state and the ground state, and the energy difference of the second and first excited 
states (Ei — E ) / (E 2 — Ei), as a function of A for Ui = 3U . The inset in Fig. 2(b) shows the 
magnification of the transition region. The figures clearly show that as A becomes smaller 
than 1, the energies of the first excited and the ground state merge together. These two 
states become quasi-degenerate, whereas the energy gap to the second excited state remains 
finite. Since the ground state is even, and the first excited state is odd with respect to the 
A — B symmetry, and since they both are Schròdinger cat states, their sum or difference 
describe the "dead" or "alive cat", respectively. 

In Fig. 3(a) and 3(b) we have plotted the energy of the first, second and third excited 
states with respect to the energy of the ground state as a function of A for A^ = 1000 and 
N = 10000, respectively. This figures clearly illustrate that, as expected, merging of the 
energy levels occurs not only for the two lowest ones, but also within consecutive pairs of 
levels, i.e. E 3 becomes practically equal to E 2 for A < 1, etc. 

Finally, in Fig. 4 we have plotted the A-atom number distributions for the ground state 
(i.e. the coefficients {q^l 2 front Eq. (|23"D as a function of m) for N = 1000 [Fig. 4(a)] and 
N = 10000 [Fig. 4(b)] for different vàlues of A. These vàlues belong to the transition regions 
between the dashed lines in Figs. 3(a) and 3(b). 

The comparison of these results with the mean field theory and its improved version 
is very satisfactory. Mean field theory is practically exact outside of the transition region, 
and gives errors of the order of 0(1/N). The improved mean field approximation of the 
subsection III. B does the similar job for all vàlues of the parameters, i.e. including the 
transition region. This result indicates that a similar improved mean field approach can be 
used for the complete field theoretical model. We adopt this approach in the next section. 

Finally, the results indicate that due the finite energy gap between the ground and first 
excited state it is possible to prepare and detect the Schròdinger cat state in the following 
manner. Obviously, direct cooling of the system to the absolute ground state, which for 
A < 1 is a cat state, would be a hardly possible task. The idea is therefore to first cool the 
system to a temperature T close to zero (i.e. such that (Ei — E )/k B T < 1) for A > 1. Note, 
that this is only possible in this regime of A, since only there the first excited state energy 
is high enough, so that practically all of the àtoms can be cooled down to the ground state. 
Then we decrease A (À) adiabatically and enter the "Schròdinger cat phase": the system 
remains in the ground state, which now becomes the "Schròdinger cat" state. 

Infernal state selective atom counting would then reveal a two-peaked structure corre- 
sponding to the \ipfi) component. This, of course, would not prové the coherence. In the 
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most general case, this would require tomographic techniques to reconstruct the complete 
density matrix of the two mode system, similar to those developed for photon fields [30 . 



IV. ANALYSIS OF THE QUANTUM FIELD THEORETICAL MODEL 

Here we analyze the full problem described by the Hamiltonians ([!]) and (|j) that ac- 
count for the atòmic motional degrees of freedom. Given the similarities of this model to 
the two-mode model analyzed above, we follow similar steps as in section III. In the first 
subsection, we apply the mean field approximation to characterize the ground state of the 
full Hamiltonian. In principle, the exact solution of the mean field equations is already 
very difficult to handle and requires numerical treatment. We have used instead two differ- 
ent methods to analyze it: the Thomas-Fermi approximation in subsection IV.C, and the 
Gaussian variational ansatz for the single particle wavefunction in subsection IV.D (for both 
methods c.f. |]2H||). In both cases we find qualitatively the same results as for the two-mode 
model; in particular, in the strong interaction and low intensity limit (the case b.2 above) 
there are two degenerate solutions of the mean field equations. In the subsection IV. E we go 
beyond the mean field approximation to analyze these results. Finally, in the last subsection 
we utilize an even more sophisticated model to approximate numerically the eigenstates of 
the system. 



A. Mean Field Approximation 



As in subsection III. A, we assume that the ground state of the system is a state for which 
all the àtoms are in the same single particle state 



(r\i) 1 )=a 1 {r)\A)+{3 l {r)\B), 



where 



i. 



(24) 



(25) 



The collective ground state of the whole system will then be, using the second quantization 
description, 



M = |V>i> ® IV>i> • • • l^i) 



7 iV! 



d 3 r 



N 



|0>, (26) 



where |0) denotes the vacuum state. The mean energy of this state can be easily calculated, 



E(a,a*,[3,P*) = ^N\H\ÏJ N ) 



d 3 ra(r)' 



V 2 r 2 U ( 
"T + 2 + ^ 



Ul 



d 3 r[a(r)*p(r){3(r)*a(r)] + (a «-> (3). 



a(f) 
(27) 



Here, as in the case of the two-mode system, we have defined Üç, = Uq(N — 1)/N, Üi = 
Ui(N — 1)/N (for simplicity we will omit the tilde in the following), and a(r) = \/Nai(r) 
and f3(r) — \f~Nf3\{f). The normalization condition (E5f) becomes now 
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I d 3 f(\a(r)\ 2 + \(3(r)\ 2 )=N. 



(28) 



In Eq. fl2~7|) the expectation value of the energy is expressed as a functional of the single 
particle wavefunctions a(r) and (3{r). The goal is now to minimize this energy with respect 
to these functions. In general, it is a difficult task that can be treated only numerically. 
In the next subsections we will follow two different approaches to find the solutions to this 
problem: first, we will analyze the Thomas-Fermi limit, and second we will use Gaussian 
ansatz. 



B. Thomas Fermi approximation 



In order to minimize ( p?| ) we calculate the functional derivatives of the mean energy 
E(a, a*, (3, (3*) with respect to a, (3, and their complex conjugates using a Lagrange multiplier 
fi to ensure that the normalization condition fl2"8|) is fulfilled. This leads to a set of coupled 
nonlinear Schodinger equations, 



~+Ç+uoHi f )\ 2 +u 1 \m 2 

^ + Ç + U \(3(r-)\ 2 + UMrl\ 2 



a(r) — \(3{r) = /xa(r), 
(3{r) - Xa(r) = fi(3{r). 



(29a) 
(29b) 



These equations are equivalent to Eqs. (]T0|) for the two-level model. In the Thomas-Fermi 
approximation one assumes that the terms involving V 2 can be neglected in comparison 
with the interaction and potential terms. 

According to Eqs. (|2"9"|), for À ^ 0, at any position r, if (3{f) = then a(f) = 0. This can 
be understood as follows: consider, for example, that at some point a(r) ^ and j3(f) = 0, 
then the làser will take partides from the state \A) to the state \B), which will imply that 
(3{f) 7^ 0. This is not the case if À = 0, where displaced solutions in the Thomas-Fermi limit 
are indeed possible f2(|. Thus, we can concentrate on the positions r where a(r),j3(r) ^ 0. 
Dividing (p9|) by a(f) and (3(r) respectively, and taking the difference we find 



U -U! + 



X 



a{r)(3{r) 



\a(f)\ 2 - (3(f)\' 



0. 



(30) 



which resembles very much Eq. (|TTJ) . The analysis is, however, a bit more complicated. As 
before, there are two kinds of solutions, \at(f)\ = \/3(f)\ and \cx(r)\ ^ \(3(r)\, where the latter 
exists for sufficiently small À only. In more detail, we can distinguish the following cases: 

(a) Weak interactions between àtoms in the states \A) and \B): In this situation Uq > U\ 



and the mean field wavefunction for the ground state is (|26|) with 



a(f) = (3{r) 



2(U + U x ) 



(31) 



where (for an isotropic trap in 3D) 
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r 



^N(U + U X ) 



1/5 



(32) 



(b) Strong interactions between àtoms in the states \ A) and \B): In this situation U\ > Uq, 
and one has to distinguish two cases: 



(b.l) Strong làser case: For 

2A 



A 



E/i-E/o 



> Ar 



15N 



in 



2/5 



(u 1 + u Q )-^ 



(33) 



the mean field ground state wavefunction is the same as in the case (a). 

(b.2) Weak làser case: For A < A there exist two degenerate solutions \tp N ) for 
the mean field ground state wavefunction of the form (|2ü|) with the coefficients 
a(f), and (3(r) given by 



a±(r) 



1 

V2 



A 2 



,1/2 



(34a) 



for r < ri = J2(fi — AUq), and 



au. r 



a_ r 



Ui-Uo 



Wif) + A/2] , 



(34b) 



for r 2 = y2/x + (C/i — C/o)A) > r > ri, Here, p At (r) = (// — r 2 /2)/f/ , and the value 
of the Lagrange multiplier /x has to be found by imposing the constraint (p8|) . 

Apart from numerical factors and different scalings, the Thomas-Fermi approximation 
çives results which are qualitatively similar to those found for the simple two-level model. 



C. Gaussian Variational ansatz 

The Thomas-Fermi solution found in the previous Section is vàlid for — > oo (or strictly 
speaking NUo t i — > oo), and predicts the existence of degenerate solutions for sufficiently low 
làser intensities. It is interesting to see if the effects remain for finite N. This can be 
analyzed using a Gaussian variational ansatz for the wavefunctions, i.e. setting 

a(r) = v / Ar r2/(4a) , (35a) 
/3(f) = VBe- r2/{4b \ (35b) 

with the variational parameters A, B, a, b. These parameters are not completely indepen- 



dent, since the normalization (|2~8 ) requires 



N = (2nf /2 (Aa 3 + Bb 3 ). (36) 
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Substituting this ansatz in Eq. fl27|), we find that the expectation value of the energy 
depends on the variational parameters a, A, b, B. We minimize it then with respect to those 
parameters taking into account the normalization condition (pü[). 

We have not found analytical solutions in this case. However, we have solved the problem 
numerically, and found the same qualitative results as in the Thomas-Fermi approximation. 
That is, in the case of weak interactions between àtoms in the states \A) and \B) (Uq > U±) 
there exists only one solution which corresponds to A — B and a = b. In the case of 
strong interactions between àtoms in the states \A) and \B) (Ui > Uq), for a given number 
of partides N we find that there exists A (N) such that if A > A the minimum energy 
corresponds to the case A = B and a = b again; Conversely, if A < A there exist two 
solutions with A ^ B and a ^ b. 



D. Beyond the Mean Field 



For our choice of parameters, the Hamiltonian ([!]), with @ is invariant under the op- 
eration Tab that exchanges the infernal level \A) with \B). Thus, the ground state of the 
Hamiltonian has to be an eigenstate of this operator. If we impose this condition on the 
ansatz fl2E| ) we find that \oc(f)\ = \/3(f)\. As we have seen in the previous subsections (case 
b.2), for a given number of partides there exists a certain A such that if A < A there are 
two wavefunctions of the form (^) \ip±) with a±(r) ^ /3±(r) that have lower energy than 
that given by the solution \a(r)\ = \/3(r)\. This implies in turn that there exists a better 
variational ansatz to the problem, namely 



where 



1 



N 



(37) 



(38) 



The corresponding energy is given by 

{r N \H\^ N ) ± ^ + N \H\r N ) 



l±(^l^) 



(39) 



It can be easily checked that E + < {ip^\H\ip^) < E_. Thus, similar ly to the two-level 
model, the proper ground state ansatz has the form of a Schròdinger cat state. 



E. Approximate Numerical Solution 

It is possible to use an even more general ansatz to generate better approximations to 
the real ground state of the Hamiltonian ([]]). In that case no analytical approximation is 
possible, and one has to restrict oneself to numerical evaluations. In any case, one can check 
whether the existence of Schròdinger cat states is compatible with these more elaborated 
calculations, and on may confirm the mean field solutions. In the following we use the ansatz 
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N 

l*> = E 



o Jm\(N — m)\ 



d 3 ra m (r)Í> a (r) ] 



d 3 r(3_ 



N—m 



N-m 



|0>, (40) 



where g m 's and the wave functions a m (f), and j3 m {r) are the variational parameters. To 
conform to the symmetry of the full Hamiltonian we impose additionally that 



pm(r) = a m {r), 



(41) 



If the ansatz (fï(]) is used to calculate the expectation value of the Hamiltonian (^\H\ty), 
one finds a rather complicated expression involving the expansion coeficients q m and the 
wavefunctions a m (r). an infinite set of nonlinear Schròdinger equations that couples q m — > 
Qm, 9m±i as we U as a m «mi «mil, ot-N-mi a N-m±i- A solution of these equations seems to 
be an impossible task, but fortunately the equations simplify in the limit of sufficiently large 
N. We have proved using the systematic \/y~N expansion, that in this limit one can simply 
substitute m ~ m ± 1, which implies a m ~ a m ±i, as well as q m ~ q m ±i in the equations for 
Oí m s. 

The resulting set of differential equations for a m (r)'s has the following form: 



V 2 r 2 



— + — + U \à m (r)\ +U 1 \à N - m (r) 



a m (r) - 2Xa N - m (r) = /j, m a m (r) 



(42) 



(43) 



is fulfilled. The above equations have to be accompanied by the linear Schròdinger equations 
for q m that contain tri-diagonal coupling to q m ±i- The coefficients in the latter equations 
depend, however, functionally in a highly nonlinear way on the a m 's, 



where ot m {r) = ^Jma m {f) , and with /i m such that the normalization condition 

m = d 3 f\à m (f) | 2 



Eq m — E m q m — XL m q m ^i — XK m q m+ i, 
where E denotes the eigenvalue we search for, whereas 



E m = \j d 3 ra m (r)* 



Ui 



Ui 



tt + ~W + -rr\ a m(r)\ + ^r\ a N-m(r)\ 



^Jm(N - m + 1) / d 3 ra m _i(f)*aAT_ m+ i(f), 



y/(m + l)(N -m) / d 3 fa m+1 (r)*a N ^ m ^ 1 (r). 



The coefficients q m are normalized as 



N 



1. 



(44) 

a w (r) + (m <-> — m), (45) 

(46) 
(47) 

(48) 



m=0 



Note that Eqs. fl42 ) are similar to Eqs. (p9|) except that they take fully into account 
the change of the form of the wave function with m. Unfortunately, these equations are still 
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very difíicult to solve, even using the Thomas-Fermi approximation. We can show, however, 
that the condition ajsi- m (r) > a m (f) (m < N/2) is fulfilled everywhere, a conclusion that 
could also be reached in the context of our mean field theory. Guided by this observation, 
we have used simple Gaussian functions to approximate the solutions of Eqs. (|^). In other 
words, we have set 

a ro (f) = {A~ m e-^l^-\ (49) 

and minimize the mean field energy numerically with respect to the a m . Note, that the 
normalization condition implies automatically that m = A m (27ca^ n ) 3 ^ 2 , so that the value of 
A m is determined by the value of a m . In an even more sophisticated attempt we have used as 
a variational ansatz a sum of two different Gaussians of the form (f4*D"|). This calculation has 
led to practically the same results as the ones obtained with a single Gaussian ansatz. For 
this reason we present below numerical results corresponding to a single Gaussian ansatz. 

Our main results are shown in Figs. 5 and 6. Fig. 5 is the straightforward analog of 
Fig. 2(a). We have plotted there the ratio between the energy difference of the first excited 
state and the ground state, and the energy difference of the second and first excited states 
(Ei - E )/(E 2 - Ei), as a function of A for iV = 1000 and \J X = 3U . We have used the 
parameters: a s ^ = a# = 50 nm, a s ^ B = 150 nm, xq = 3 /im, Two = 100 Hz. The figure clearly 
shows that as in the the case of the two- mode model as A becomes smaller than 1, the 
energies of the first excited and the ground state merge together. These two states become 
quasi-degenerate, whereas the energy gap to the second excited state remains finite. 

Similarly, Fig. 6 is an analog of Fig. 3(a). There we have plotted the energy of the 
first, second and third excited states with respect to the energy of the ground state as a 
function of A. As already seen in the two-mode model, merging of the energy levels occurs 
also among consecutive pairs of levels, i.e. E 3 becomes practically equal to E 2 for A < 1. 

In conclusion we stress that the similarity between the Figs. 2(a), 3(a) and 5, 6 is 
remarkable. Clearly, the complete field theoretical model leads to the same physics as the 
two-mode model. As A is adiabatically decreased, the system enters the "Schròdinger cat 
phase" in which the ground state is a linear superposition of two states for which the A-atom 
number distributions are significantly distinct. 



V. IS A SCHRÒDINGER CAT STATE EXPERIMENTALLY FEASIBLE? 

We conclude with a summary of requirements to observe Schròdinger cats in an exper- 
iment. While these necessary conditions to prepare and preserve cat-like states are not 
fulfilled in the present generation of Bose-Einstein experiments |23]], they might provide a 
guideline and motivation for future experimental work. 

(i) A — B symmetry. The calculations in the present paper assume an A — B symmetry. 
We believe that this assumption is mainly a technical point in the theoretical calculation, 
but discuss this now in more detail. 

The choice of equal scattering length, reasonable, and agrees with the 



recent theoretical calculations P5| . The assumption of equal trap freqüències, uja = ^b, 
however, is typically not fulfilled in magnètic traps, since àtoms in different infernal states 
feel different (magnètic) potentials, but could be realized in principle in an optical dipole 
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trap, where it is assumed that the two states have the same electron configurat ion, so that 
the far off resonant làsers induce the same lightshifts. Another way of achieving equal trap 
freqüències, is to compensate an asymmetry uja ^ ojb using an appropriately detuned làser. 
Such compensation is exact in the case of the two-mode model. In the case of the complete 
field theoretical model it requires a little more care. The idea is that the necessary and 
suficient condition for a ground state of the system to be a Schròdinger cat state, is that 
there exist two distinct and degenerate minima of the energy in the mean field approximation. 
If uja 7^ ojb the mean field theory would typically lead to two minima of the energy function, 
but with slightly different energies. The reader can easily convince him/herself that this is 
the case for the two-mode model, while the analysis for the complete model is more technical, 
but otherwise analogous (especially in the Thomas- Fermi limit). That means that in such 
cases we would have one global, and one local minimum of the energy, and neither of them 
would correspond to a "Schròdinger cat" state. They will, however, be characterized by 
different numbers of àtoms in the states \A) and \B). The point is that the energies of these 
minima, and the corresponding A- and 5-atom numbers can be deformed in a continuous 
manner by changing A. At so me point one arrives at the situation when the two minima 
become degenerate, and at which the true ground state becomes a linear combination of the 
two, i.e. becomes a "Schròdinger cat" . 

(ii) Instability condition a s ^ B > a s ^ B . The existence of cat states imposes the (instability) 
condition a s ^ B > a s ^ B (i.e. U\ > U ). As we know, the present experiments |23| with Rb 
àtoms allow for simultaneous evaporative and sympathetic cooling because the inelastic 



collision rates are small, which is directly related to the fact that a s X B — o s a — Qr P6 
While the condition seems not to be satisfied for atom species used in present magnètic 
trap experiments, we stress that future experiments might be based on different traps, for 
example, far-off resonance traps using hightly detuned làser light. This will open up the 
possibility of trapping new infernal atòmic states. Consider for example a total angular 
momentum F — 1, and condensates in the states mj? = ±1. Furthermore, we assume 
that the level mj? = does not participate in the collision dynamics (this can be done, 
for example, by shifting it with a làser). If the singlet scattering length is larger than the 
triplet scattering length then the condition a s ^ B > a s ^ B will be fulfilled |[28|| . In addition, in 
principle, to modify the atom-atom scattering lengths [ 27fl . 

(iii) Cooling to the ground state. Preparation of a "cat" requires the cooling to the 
ground state of our system, i.e the preparation of a pure state. The sufficient condition is 
that the temperature has to be k b T < Ei, where E% is the energy of the first excited state 
of the total Hamiltonian (for example, in the ideal case, that would require that more or 
less N — 1 partides are in the ground state, and one is in the first excited state). We stress 
that this requirement is much stronger than the requirement of having most of the àtoms in 
the single particle ground state, i.e obtaining a macroscopic occupation of the ground state, 
as observed in current BEC experiments ||T] [| . 

We illustrate this by an example: if we have N partides and 80% of them are in the 
ground state, the population of the collective ground state is only 0.8^. Thus, the tem- 
peratures required to observe a "cat" are much lower, such that practically all àtoms are 
in the collective ground state. Wether existing cooling techniques, in particular evapora- 
tive cooling, can be extended into this regime remains to be investigated. However, the 
fact that we are dealing with bosons instead of distinguishible partides helps significantly 



16 



to prepare a pure "cat" ground state. Consider a one dimensional situation. For N dis- 
tinguishible partides, there are N possible collective excited states with the same energy. 
Thus, the temperature required to have most of the population in the collective ground state 
is kbT < Ei/N. On the contrary, for bosons this temperature is N times higher, kT < E\. 

(iv) Decoherence. Finally, we should address the question of decoherence. Obviously, 
atom losses (such as those due to inelastic col·lisions) would destroy the "Schròdinger cat" 
state very rapidly. In fact, in the extremal case when one has the cat |0, AT) + \N, 0), 
already one atom loss would be enough to distort completely the coherent superposition (c.f. 
|T^,TPJT^]). We stress, however, that the situation here is similar to that of the experiment 
of Brune et al. ||15|| . The "cats" that live long enough to be observed must be mesoscopic. 
In fact, the "Schròdinger cat" states displayed in Fig. 4 belong to that category. They allow 
for loss of many àtoms without the complete smearing out of their quantum mechanical 
coherence. If they are created, they could allow for the study of the gradual decoherence 
process, in a similar manner as has been done in Ref. fTSfl . 
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FIGURES 




FIG 1. The ground state energy Eq as a function of A for N = 1000 and U\ = 3Uo calculated 
for the two-mode model. 




A/N 

FIG. 2. (a) The ratio between the energy difference of the first excited state and the ground 
state, and the energy difference of the second and first excited states {E\ — Eq)/(E2 — Ei), as a 
function of A for U x = 3U and for N = 1000; (b) Same as (a), but for N = 10000. The inset 
shows the magnification of the transition region. 
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0.8 0.9 1 1.1 



À/N 

FIG. 3. (a) The energy of the first, second and third excited states with respect to the energy 
of the ground state as a function of A for N = 1000, and other parameters the same as in Fig. 1. 
The dashed lines mark approximately the transition region; (b) Same as (a), but for N = 10000. 
The inset shows the magnification of the transition region. 
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500 1000 5000 10000 

FIG. 4. (a) A-atom number distributions for the ground state (i.e. the coefficients \q%\ 2 from 

Eq. ( |23| ) as a function of m) for A" = 1000 for the vàlues of A indicated. These vàlues belong to 

the transition region between the dashed lines in Fig. 3(a); (b) Same as (a), but for N = 10000. 

The vàlues of A belong to the transition region between the dashed lines in Fig. 3(b). 








0.6 0.7 

FIG. 5. The energy difference of the first excited state and the ground state, and the energy 
difference ofthe second and first excited states {Ei—Eq)/{E2 — E{), as a function of A for N = 1000 
and U\ = SUo, calculated for the complete quantum field theoretical model. The parameters are: 
a A =a B = ^0 nm ) &AB = 150 nm, x = 3 jttm, fiu = 100 Hz 
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0.6 0.7 

FIG. 6. The energy of the first, second and third excited states with respect to the energy of 
the ground state as a function of A calculated for the complete quantum field theoretical model. 
The parameters are the same as in Fig. 5. 
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